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The f{R) gravity models proposed by Hu-Sawicki and Starobinsky are generic for local gravity 
constraints to be evaded. The large deviations from these models either result into violation of local 
gravity constraints or the modifications are not distinguishable from cosmological constant. The 
curvature singularity in these models is generic but can be avoided provided that proper fine tuning 
is imposed on the evolution of scalaron in the high curvature regime. In principle, the problem can 
be circumvented by incorporating quadratic curvature correction in the Lagrangian though it might 
be quite challenging to probe the relevant region numerically. 



PACS numbers: 98.80 Cq 

I. INTRODUCTION 

The growing faith in the late time cosmic acceleration is directly supported by observations of high red-shift 
supernovae and indirectly by observations on microwave background, large scale structure and weak lensing. What 
causes the repulsive effect, in the cosmic expansion, is one of mysteries of modern cosmology at present. Theoretically, 
the phenomenon can be accounted for either by supplementing the energy momentum tensor by an exotic matter 
component with large negative pressure {dark energj/) P, [3, H, UH, 0, 0] or by modifying gravity itself. Cosmological 
constant, the simplest candidate of dark energy, is plagued with fine tuning problem of an unacceptable level Q. Scalar 
field could provide an interesting alternative to cosmological constant Q. They can mimic cosmological constant like 
behavior at late times and can give rise to a viable cosmological dynamics at early epochs. Scalar field models with 
generic features are capable of alleviating the fine tuning and coincidence problems [3|. As for the observations, at 
present, they are absolutely consistent with A but at the same time, a large number of scalar field models are also 
permitted. Future data should allow to narrow down the class of permissible models of dark energy. 

As an alternative to dark energy, the large scale modifications of gravity could account for the current acceleration 
of universe. We know that gravity is modified at short distance and there is no guarantee that it would not suffer 
any correction at large scales where it is never verified directly. Large scale modifications might arise from extra 
dimensional effects or can be inspired by fundamental theories of high energy physics. On purely phenomenological 
grounds, one could seek a modification of Einstein gravity by replacing the Ricci scalar in Einstein-Hilbert action by 
However, any large scale modification of gravity should reconcile with Local Gravity Constraints and should 
have potential of being distinguished from cosmological constant. Since the general theory of relativity is in excellent 
agreement with local gravity phenomenon, it is quite challenging to construct a viable model of f{R) gravity along the 
said lines. Stability requires that the first and the second derivatives of /(i?) with respect to the Ricci scalar R should 
be positive definite. Most of the corresponding modifications of the Einstein-Hilbert action are either cosmologically 
un- viable or can not be distinguished for the cosmological constant. 

The class of models proposed by Hu-Sawicki and Starobinsky (HSS) is of great interest. These models can evade local 
gravity constraints and have potential capability of being distinguished from the cosmological constant [HI (see also 
Ref . j 12| ) . However, they are quite delicate — the minimum of the scalaron (scalar degree of freedom present in f{R) 
gravity) potential which corresponds to dark energy in these models is very near to field configuration corresponding 
to infinitely large value of R for solar physics constraints to evaded. Thus it is quite likely that the scalar field, which 
controls the space-time curvature, hits singularity while evolving near the de-Sitter minimum 

The problem becomes acute in high curvature regime but can be circumvented by carefully tuning the parameters of 
the modelfli,!!^. 

The HSS models are characterized by a finite potential barrier between the minimum of the scalaron potential and 
the curvature singularity and hence are vulnerable to singularity. Recently, the HSS models were modified such that 



the said potential barrier is infinite and the curvature singularity is hidden behind the infinite potential barrier [21 1. 

In this paper, we examine the deformations of HSS models and demonstrate that these models are generic to local 
gravity constraints. We also argue that the viable resolution of curvature singularity can be provided by adding higher 
curvature terms to the originally proposed form of f{R) in Refs.p^. 
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II. LARGE CURVATURE SINGULARITY VERSUS THE LOCAL GRAVITY CONSTRAINTS 



The action of f{R) gravity in Jordan frame in the presence of matter described by the matter Lagrangian Cm is 
given by 



S 



fiR) 



(1) 



where the matter Lagrangian depends on the metric g^j, and the matter fields. In what follows, it would be convenient 
to write /(i?) in the following form, 



/(i?) = i? + A, 



(2) 



where A describes the correction to Einstein-Hilber action and A^^ denotes its derivative with respect to the Ricci 
scalar R. The f{R) theory apart from the spin-2 object necessarily contains a scalar degree of freedom which becomes 
clear either by taking the trace of the modified Einstein equations obtained from ([1]) or by passing to the Einstein 
frame. Indeed one can always make a conformal transformation which converts the original action ([T]) into Einstein- 
Hilbert action along with a canonical scalar field (j) which directly couples to matter. The solar system and equivalence 
principle bounds give a strong constraint on the magnitude of the scalar field in the Einstein frame. The potential 
of field (j) is uniquely constructed from the Ricci scalar R. 

We now transform the metric using the conformal transformation. 



The action in the Einstein frame is given by 



where the coupling gc and V are given by, 



9c 



1 

V6' 



V 



Rf.R - f 



As shown in Ref . |22l| . the thin shell parameter is given by 

Afe _ (j)B ~ i'A 
fc &9c^c 

where <j)A, 4>b are corresponding to the minimum of the effective potential 

K//(</') = ^(0) + eS=V, 



(3) 



(4) 



(5) 



(6) 



(7) 



inside and outside the spherical body respectively and f&c is the gravitational potential of the test body (Sun/Earth). 
Let us consider the variants of HSS models plj|. 



a(3Rc 



-1/13 



Rc > 



(8) 



The conditions for the cosmological viability of f{R) models can be understood by considering two quantities [22 



Rf,RR Rf.R /f-s 

m=— , — (9) 

J.R J 



The presence of a viable saddle matter era demands that 



m(r «-l)wO, m'(rw -!)>-! (10) 



3 



The conditions (fTO]) are satisfied for the model ^ provided that 



n > 0, and (/3 > 0, or /3 < ~n) (11) 

In what follows we shall consider the case of (n, /3) > 0. In fact we know [l^l that for n > and (3 < the model 
is not distinguishable from the ACDM. 

Let us emphasize that HSS models in Starobinsky parametrization corresponding to n = 2 and (3 < 1 has a moderate 
dependence on R allowing the local gravity constraints to be evaded. 

Let us now analyse extended HSS niodels[2l[ described by ([8]). In this case, in the high curvature regime R ^ Rc, 
we obtain, 

A.fl «-an(^|-) ^ \ (12) 

which shows that A ^ 1 in the case under consideration (n, (3 > 0) for moderate values of a. Using expression for 
(j) given by Eq.® and the fact that Aji ^ 1, wc find, 

<^=y|ln(l + A,«)«^A,^, (13) 
Wc next estimate R corresponding to minimum of the effective potential. 



R{1-A^r) + 2A 



9ce 



9^4' p* = (14) 



(1 + A.i?.)2 

which simplifies in case of the generic approximation, A^^ 1, A ^ i? and gives rise to following expression for 0„im 

0mm ~ —A^nln^p, « — —an — . (15) 



2 2 \p* 

Hereafter, we shall use the notation p for matter density instead of p* in Einstein frame. From the fact that pA, 
the energy density inside the test bodies (Sun/Earth) is of the order of 1 g/cm'^ which is much larger than the 
density outside {pB ~ 10^^^ g/cm^ of the baryonic/dark matter density in our galaxy), it follows from Eg. pS]) that 

\M < |0b|, 

(PB \PaJ 

which allows us to write the thin shell condition in the convenient form 

|0b| =^ V6$,^, (17) 



< 



5.97 X 10~" (Solar system test), 
3.43 X 10"^^ (Equivalence Principle (EP) test) 



(18) 



We have used ^ < 1.15 x 10"^ $c ^ 2.12 x 10"^ for the Sun and < 2 x 10"^ <i>c ~ 7 x 10"^" to respect the 
equivalence principle constraint. In what follows, we shall investigate the modified HSS models ([8]) for different values 
of model parameters. Let us first consider the case oi P oo and n ~ I 2l| . 

The de-Sitter minimum in free space is given by 
which gives rise to the following relation for a 



a;i = ^ (21) 



-xi +2(l + a;i)ln(l + .Ti)' R^ 



Relation (j21[) implies that a is always positive definite for any xi and that a ^ 1 as cci ^ 0. For moderate values of 
a, the de-Sitter minimum corresponds to Ri ^ pc {pc — 10~^^ g/cm~^). For instance, in case of a = 2, we find that, 
Ri ~ 6i?c which gives the estimate for \(f)B\ as |0b| > 10~^. This is clearly ruled out by the thin shell condition 
We next investigate the model ^ for arbitrary values of parameters. 
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A. Constraint for the general /?, n, a 

Let us define the dimensionless variable x as x = R/ Rc and write the expression of interest in terms of x, 



A 



-a/3i?,{l-(l + x")-i/^}, (22) 



A,R = -naa:""i(l + a;")-i/'^-\ (23) 
_ RA,n - A 
- 2(1 + A,«)2' 

yRc (1 + x")-l//3-l _ (1 + 2;") [_1 + (1 + a;»)l//3] ^} 



2 [-1 + ?i.Q:a:;"-i(l + a;")-V/3-i]^ 

The dc-Sittcr minimum in free space in this case corresponds to 

xi(l + x5')i+i/'3 



(25) 



(26) 



-nx'l + 2{l + x'^) [-l + (l + a;^)i//3]/3' 
For /3 ^ 00, these equations reduce to 

A = -ai?cln(l + a;"), (27) 



\r = — -, (28) 

x" + 1 

V = «i?c x^jl + x") [nx" - (1 + rr") ln(l + x")] ^ 



Xl(l+X?) 



^?i.T'i' + 2(1 + x'l) ln(l + x^*) 



B. case: /3 — > cxd and n > 2 



(30) 



Our numerical analysis shows that a is positive definite for all values of xi provided that n < 2. However, for larger 
values of n there exist values of xi for which a is positive. In case of n ^ 10 corresponding to xi > y/e = 1.649, the 
parameter, a is always positive as shown in FigUJ In this case x" ^ 1 and we obtain 



xi(l + x5') x"+^ Xl 



-nx;^ + 2(l + .x'/)ln(l + x'i') x']'(21n(x'i') -n) ?i(21nxi-l)^ 
which we shall use to confront the model with solar tests 



(31) 



\^B\^~^AM,^,^=^jP^-^na^^J -^na(x,fj - ^- x 10^ (32) 



I na 



Substituting, a from Eg. pTj) in (|32p . we have 

\M « ^ . X 10-5 (33) 
2 mxi — 1 

To satisfy, < 0(10"^°) (EP constraint), we need Inxi to be very large number which implies that the model can 
not be distinguished from ACDM. 

C. case: /3^oo,7i^0 



xi(l + x?) _ xi(l + l) _ Xl ^g^^ 



~nx?;' + 2(l + x'/)ln(l + x5') (1 + I)ln2 ln2 
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FIG. 1; Plot of the parameter space for Ri/Rc varying from 0.01 to 100 and n ranging from 0.005 to 10. The plot shows that 
there is no region in this parameter space for the case (3 ^ oo for the local gravity constraints to be satisfied, \4>b \ ^ 10"^". 



and 



I0f 



%/6 ( Pb 



\/6 na 



X 10" 



-\/6 n 



X 10 



-5 



(35) 



then n < 0(10 ^''), as implied by the EP constraint, which makes the model indistinguishable from cosmological 
constant. 



D. 



case: /3 — > cxD,n < 2 



As shown above, for arbitrary value of n, the parameter a is positive definite provided that n < 2 for all values of 
xi. This can easily be demonstrated analytically in the limit of xi 0, 



Xi 



„1-1 



-nx'1 + 2(1 + .t'/) ln(l + x'l) -~nx'1 + 2{l){x'l) 



(36) 



In the region of positive a and 0.2 < n < 2, our mmrcrical estimates show (see FigH]) that \(f)B\ > 10 ^. As mentioned 
before, the model may be compatible with solar test for n < 10^^" but reduces to KCDM . Let us note that the class 
of models il 



-aRr. 1 



R_ 



(37) 



is practically not distinguishable from cosmological constant as the local gravity constraints impose severe restriction 
on n, namely, n < 10"^*^. 

So far we have focusscd on large P limit of Starobinsky model as singularity is clearly avoided in this case. A 
comment on the finite /? behavior of the model is in order. In this case, the analysis requires numerical treatment. 
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FIG. 2: Plot of the potential V versus tp for different values of parameters. The plot sho-ws that the potential barrier becomes 
large for the large values of /3. The de-Sitter minimum is also seen to shifts towards the singularity as /3 decreases. 



In Fig. 3, we have displayed the parameter region consistent with local gravity constraints in case of finite values of 
n and p. In agreement with Ref.jl^l, we find that the local gravity constraints are satisfied provided that n/P ^ 2, 
see Fig.3. Furthermore, in f{R) gravity, the power spectrum acquires an additional slope [ll| which is constrained in 
Ref. [25. As demonstrated by Starobinsky , n/(3 should satisfy the constraint, ^ > 4. Thus the model proposed in 
Ref . [2]| . with /3 — > 00, violates this constraint too. If we adhere to observational constraints imposed by local gravity 
constraints, the model is vulnerable to curvature singularity. It is really interesting that the height of the barrier 
between de-Sitter minimum and curvature singularity turns out to be proportional to f3 which is heavily constrained 
by local gravity constraints. In what follows, we shall address this issue. 



III. FINITE TIME GENERIC SINGULARITY AND ITS RECONCILIATION 

Let us note that in the limit of i? — !■ 00, A^^ and the maximum of the potential is located at "0 = 1 whose 
magnitude is given by 



V = 



RA^R - A 



2il + A.R) 



R=oo 



(38) 



Since limn^oa RA^r = and limn^ac A = —a(3Rc, we find that lim/j^oo V/Rc = The minimum of the potential 
in free space given by ([38)) can be estimated numerically, Vmin/Rc — C>{1) at i? ~ In this case the height of the 
potential barrier for large value of f3 is approximately equal to /?/2 as shown in the Figl^l The local gravity constraints 
impose a restriction on the height of the barrier or equivalently, the parameter /3 for a given value of n and Re- In 
case of 71 = 2, Ri/Rc = 4 and (3=1, the model passes both the local gravity constraints. For large values of /3, the 
height of the potential barrier becomes large thereby hiding the singularity but resulting into clear violation of local 
gravity constraints. We also observe that taking small values of f3, the de-Sitter minimum shifts towards singularity, 
see Figl21 This implies that we should have moderate values of parameters for a viable evolution. Situation gets 
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FIG. 3: The figure shows the allowed regions of the parameter space which satisfy the thin shell condition corresponding to EP 
constraint for the various values of /3 in case Ri/Rc and n range from 0.1 to 10 and from 0.05 to 10 respectively. It is clearly 
seen that n/ f3 > 1.7 to satisfy local gravity constraints. 



worse when we move to high density regime whose treatment requires extreme fine tuning of initial conditions of the 
field [131. 

As we have seen that the size of |<^mm| ~ \^,r\r=p for any viable f{R) gravity and is constrained to be less than 
0(10"^°). This means that the minimum of the potential corresponding to ■0 = 1 + ^,r, is very close to ■0 = 1 even 
in the case of baryonic/dark matter density. 

It should be emphasized that in case of large curvature, the quantum effects become important leading to higher 
curvature corrections. Keeping this in mind, we can incorporate iiR^/Rc term in the model under consideration [Til 
[i5|(see Ref.flBl on the similar theme) which allows us to move the singularity away from = 1. The Big Bang 
nucleosynthesis constraint at T ^ MeV {z ^ 10^°) or R ^ lO'^'^pc tells us that the correction term should satisfy the 
following condition [27j|, 

<c R. (39) 

Rc 

If we choose Rc ^ Pc, wc find that fi <C 0(10^^"). In case of neutron star with p ^ 10*'^ pc, the parameter p is 
constrained to be /i << 10^*'^. The local gravity constraints arc satisfied in this case as 

(from /xi?V^c term) ~ A^r\r^io'^r^ 2-^10^ R^ <^ O{10-^^). (40) 

Let us note that in case we intend to describe inflation with the help of R^ terms — a la Starobinsky model, the 
numerical value of p is much smaller than the quoted value. Indeed, the mass of scalaron (i?y^/6/i^/^), if it is to 
be infiaton, should be 10~^Mp[ll[ which implies that p is much smaller than its numerical value quoted in case of 
neutron star. Such a correction does not disturb the neutron star physics and nucleosynthesis constraint but can help 
in avoiding the curvature singularity. As a result, the correction term can not contribute any effect to the local gravity 
experiments. This implies that the behavior of the model improves in the high curvature regime though it might be 
quite challenging to probe it numerically. 
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r 



FIG. 4: The (r, m) plane for the HSS model (/? = l,n = 2) with the additional term -j^R^- The dashed diagonal line is the 
critical line m — —r — 1. The numerical value of a is chosen such that the condition of the stability of the future de Sitter 
stage is satisfied. 

We can sec from Fig.(|4|) that in case the BBN condition on /i is satisfied, the model would have a standard matter 
phase. Then a small /i is not only necessary for the BBN but also for the matter phase. Furthermore because of this 
additional term the curve m(r) cross the line m = —1 — r for a finite R while R = oo for the HSS model. Thus we 
can connect the early phase of accelerated expansion to the late time acceleration of universe without a singularity of 
the curvature scalar R. 



IV. CONCLUSIONS 

In this paper we have examined the variants of HSS models described by three parameters, a, P and n[2l1|. The HSS 
scenario in Starobinsky parametrization corresponds to n = 2 and /3 ^ 1. These models can satisfy the local gravity 
constraints and have potential capability of being distinguished from cosmological constant. The de-Sitter minimum 
of the effective potential of the scalar degree of freedom is quite close to curvature singularity for moderate values 
of a and /3 in these models. For larger values of a and 1/(3, the de-Sitter minimum moves towards singularity. The 
potential barrier between the de-Sitter minimum and curvature singularity is finite in HSS models which makes them 
delicate. Thus one should carefully tune the scalaron evolution such that it docs not hit singularity while evolving in 
the neighborhood of the minimum of effective potential. For a given value of n, the height of the barrier is defined by 
the parameter (3 which is large for larger values of /3 thereby hiding the singularity behind the potential barrier [2 1|. 
However, large potential barrier between singularity and de-Sitter minimum comes into conflict with the local gravity 
constraints. 

The high curvature behavior of An is extremely crucial for local gravity constraints to be evaded. In case A^ji 
slowly as it happens in case of large /3, we can not satisfy the local gravity constraints. On the contrary, if An 
approaches zero fast, the corresponding models become more vulnerable to singularity as the minimum of the effective 
potential moves very near to singularity in the high curvature regime. In this case, the models under consideration 
can hardly be distinguished from cosmological constant. 

The HSS scenario is build very carefully such that the curvature dependence of A(i?) is just right to satisfy the local 
gravity constraints and at the same time to allow to distinguish itself from ACDM. Thus the finite time singularity in 
viable f{R) models is generic. However, the safe passage of the scalaron to the minimum of its effective potential can 
be ensured by the appropriate fine tuning of the scalaron evolution poj. The fine tuning turns ugly in case of compact 
objects like neutron stars. The introduction of higher curvature terms becomes legitimate near singularity and can in 
principle improve the behavior of the model. The resulting scenario can give rise to a viable cosmic evolution. 
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